Thinning of superfluid films below the critical point 
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Experiments on 4 He films reveal an attractive Casimir-like force at the bulk A-point, and in the 
superfluid regime. Previous work has explained the magnitude of this force at the A transition and 
deep in the superfluid region but not the substantial attractive force immediately below the A-point. 
Utilizing a simple mean-field calculation renormalized by critical fluctuations we obtain an effective 
Casimir force that is qualitatively consistent with the scaling function ■& obtained by collapse of 
experimental data. 
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Experiments by Garcia and Chan [l| have challenged 
our understanding of finite size behavior in thin films of 
4 He close to the A-point. The experiments monitor the 
thickness of a wetting film of helium suspended above a 
4 He bath. Near and below the bulk transition temper- 
ature the film becomes thinner. Since the microscopic 
intermolecular interactions are not expected to be mod- 
ified at the A-temperature, the thinning can only be due 
to collective behavior at and below the transition. In- 
deed, the reduced thickness has been interpreted as due 
to an attractive, fluctuation-induced, Casimir-like force. 
Subsequent results by Ganshin, et. al. confirm that 
this force is consistent with finite size scaling, according 
to which 0, 0], for a slab of width, L, the effective force 
per unit area f(T,L), acting to thicken or thin the film, 
has the general form 

f(T,L) = ^m 1/u )- (i) 

Here, t = (T — T c )/T c is the reduced temperature, mea- 
suring the distance from the bulk critical point T c , # 
is the dimensionless amplitude of the scaling function, 
and v is the correlation length exponent. Figure [T] shows 
the experimental data for the dependence of the function 
■d(x) on the scaling combination x = tL x l v ', as presented 
in Ref. Q (see also Fig. 3 in that article). This plot 
combines the results from 4 He films of three different ini- 
tial thicknesses of 238, 285, and 340 A. It is important 
to note that these films are in equilibrium with the bulk 
Helium liquid. 

As shown in the figure, the force amplitudes extracted 
from the data associated with the three films collapse 
onto a single curve, in clear vindication of finite size scal- 
ing. The behavior of the scaling function at and imme- 
diately above the A transition is well described by the 
two-loop renomalization group calculations of Krech and 
Dietrich [1, 0, Sj- Another feature that is successfully ex- 
plained by theory is the finite negative value of d well 




FIG. 1: The Casimir force amplitude i? for 4 He films of three 
different thickness L, as a function of the scaling variable 
x = tL 1 ^, from the data by Ganshin, et. al. @, Ej] . Despite 
the differences in film thickness, the scaled forces collapse onto 
a single curve with a minimum at x = —9.7 ± 0.8A 1 /". 



below the critical point: it was demonstrated in Ref. [9j 
that a combination of phase fluctuations within the film 
and surface deformations can account for this value. 

However, the most prominent feature in Fig. [T] is the 
deep minimum of the scaling function in the vicinity of 

x = tl}l v = -9.7 ±0.8 A 1 /". (2) 

The amplitude of the force at this minimum is roughly 
an order of magnitude larger than that of the Casimir 
forces at (x = 0) and far below (x — > — oo) the transition 
and has thus far not been reproduced by any theoretical 
calculation. In particular, the dip is considerably larger 
than the one in the simulations reported by Dantchev 
and Krech fl~oj ] , or the low temperature vortex loop cal- 
culations of Williams [H H Ell • 

In this letter, we show that the origin of the large at- 
tractive Casimir-like force measured in Refs. [3,12] is 
twofold. First, the boundary conditions that apply in the 
case of a superfluid film arc Dirichlct in that the super- 
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fluid order parameter effectively vanishes at the substrate 
and at the liquid-vapor interface. This distinguishes the 
wetting film from the periodic systems considered in Refs. 

ToL 11 . [Hj]. Second, because the film is in equilibrium 
with a reservoir of bulk fluid, it is necessary to take into 
account the free energy of the bulk in calculating the ef- 
fective force, /, acting on the film. This leads to the 
following scenario on cooling the system: Immediately 
below the A temperature, ordering in the form of a non- 
zero superfluid order parameter takes place in the bulk 
reservoir of 4 He. However, boundary effects suppress the 
appearance of superfluid ordering, even at a local level, 
in the film. To take advantage of the lower free energy 
of the superfluid state, helium atoms move from the film 
to the ordered bulk, leading to additional film thinning, 
and a corresponding increase in the attractive effective 
Casimir-like force. Upon further reduction of tempera- 
ture, local ordering finally occurs in the film; its free en- 
ergy effectively "recovers" to a value similar to the bulk, 
resulting in a rebound of the film thickness. 

Both of these contributing factors are well-displayed in 
a simple mean field treatment augmented by the renor- 
malizing effects of critical fluctuations at the onset of 
order in the bulk. Calculations based on this approach 
indicate that the minimum of the scaling function occurs 
when the dimensionless quantity y = (L/^) 1 /", where £ 
the correlation length, is at y m i n — —tt 2 = —9.87. The 
experimental value of y* in is equal to — 7.4±0.6, as is ob- 
tained by setting t = in E q-I3 with £o = 1-2 A 

14j . We note at the outset that the mean-field calcu- 
lation, details of which follow, yields a minimum am- 
plitude which is significantly greater than the observed 
value. We are optimistic that a proper treatment which 
correctly accounts for the effects of fluctuations — and for 
the fact that the superfluid order parameter should be 
modeled by a two-component vector — will lead to a sub- 
stantial reduction in the force. However, it is clear that 
any proper theoretical model will have to embody the 
elements underlying the results reported here. 

The simplest depiction of the superfluid transition in a 
thin film geometry is based on the Ginzburg-Pitaevskii 
Hamiltonian[l5|. which follows from the Ginzburg- 
Landau model[16]. In one dimension (moving through 
the film) the effective Hamiltonian is 



n 




+ -(f>{z) 2 + u4>(z) 4 > dz, (3) 



where <j)(z) is the order parameter, r is the reduced tem- 
perature, and u is the fourth order coupling constant, 
which quantifies the system's self interaction and pro- 
vides the mechanism by which ordering saturates. 

The extremum (saddle point or mean field) equation 
satisfied by the order parameter is 



= 
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d 2 (/>(z) 
dz 2 



+ rcj>(z) + Aucj)(zf 



We assume Dirichlet boundary conditions, 0(0) = 
4>(L) = 0, at both sides of the film consistent with the ex- 
perimental situation. For the bulk (L — > oo), or in a peri- 
odic and hence homogeneous system, a non-zero solution 
appears as soon as — r > 0. For the film, it is easy to see 
that a non-zero solution appears only for — r > Ao, where 
Ao = tt 2 /L 2 is the smallest (in magnitude) eigenvalue of 
the Laplacian with Dirichlet boundary conditions. 
By quadratures, the full solution to this equation is 



d(j)(z) 
dz 



(z) 2 )+ U (0(z) 4 -^), (4) 



where we have obtained Eq. ((4]) from Eq. ([3]) by exploit- 
ing the mathematical equivalence between the first ver- 
sion of the equation of state and Newton's second law for 
the motion of a particle in a one-dimensional potential. 
Equation ^ is the statement of conservation of energy 
in this context. 

In accord with the boundary conditions and obvious 
symmetry considerations, we set <fi'(L/2) — and thus 
4>{L/2) — 4> - Integrating dz/d<p in Eq. ((4]), from the edge 
of the film to its midpoint gives [TtJ 
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o A /-r(^-^)[l + ( 2u /r)(^ + ^)] 



1 K 



-r VI = V 

where K is the elliptic integral of the first kind and 
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To obtain the effective force, we calculate the deriva- 
tive of the free energy in Eq. ^ with respect to L, after 
replacing <p(z) with its extremum value. After a series 
of relatively straightforward steps, we end up with the 
following result for the L-derivative of the free energy 



dl _ 1 U<p{z) 
~dL ~ ~ 2 V dz 



z=L 



r -cb(L) 2 + uHL) 4 . (7) 



Finally, making use of the Dirichlet boundary conditions, 
and Eq. Q , we have for the derivative of the free energy 
with respect to the thickness of the system, 



dL 



■u<t>Q. 



(8) 



We note that Eq. ([5]) allows us to write 77 as a function 
of the combination y/—rL/2. Setting r = y/L 2 and mak- 
ing use of Eq. ([6]), the derivative of the free energy with 
respect to L takes the form 



~dL 



jfriiV^vMil-vW^v/Z))- (9) 
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Except for the quantity u and the prefactor of l/L 4 : 
which is appropriate to a four dimensional system (and 
is consistent with the fact that mean field theory agrees 
with 4-dimensional hyperscaling); we have a free energy 
derivative that depends on the scaling variable y. Within 
mean field theory, the quantity u is material-dependent. 
However, the existence of a stable renormalization group 
fixed point on the critical hypersurface implies a univer- 
sal value for this coefficient [1 81 ] . 

An important feature of the experiments in Refs. [HQ 
is that the wetting film is in equilibrium with the va- 
por above a bulk reservoir of 4 He. Thus, to determine 
the layer thickness we should calculate the change in free 
energy as some fluid is removed from the film and simul- 
taneously added to the bulk. The mean field calculation 
of the bulk free energy is straightforward, obtained from 
the previous result in the limit of L — > oo. The resulting 
expression for the force is 
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Figure [2] displays the above mean field scaling function, 
d(y) = fL /ksTc, as a function of y = rl? = (L/£) 2 , 
where £ is the correlation length as expressed before. Also 
shown in the figure is the data in Fig. [TJ Note that in 
Fig.[T]the scaling function is plotted vs. x which has units 
of A 1 /", while in Fig. [2] the horizontal axis represents the 
dimensionless quantity y — (L/^) 1 ^. The amplitude of 
the mean field function depends on the parameter u, and 
has been adjusted to give the observed amplitude at the 
minimum. Note that at this coarse resolution, the mean- 
field curve captures the main trends of the experimental 
data below the critical temperature. The most significant 
discrepancy is that the measured force does not go to zero 
for y — > — oo; this is due to the absence of Goldstone and 
surface modes in the mean- field analysis Q. 
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FIG. 2: The mean field scaling function t?(j/) derived from 
Eq. (|10|l (dashed curve), compared to the re-plotted experi- 
mental data from Fig. [1] The horizontal axis is dimensionless 
and provides a check of the approximations, while the vertical 
axis is adjusted so that the minima occur at the same point. 



To obtain the correct scaling of the force / with the 
film thickness L, we utilize a version of mean field the- 
ory that takes into account the renormalizing effects of 
critical fluctuations (see Ref. [Hj]). The relevant free en- 
ergy density, after renormalizing to cell blocks larger by 
a factor b = e l " , is 
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-<t>(l*f+u{l*)ci>{l*) 4 



(11) 



Close to a fixed point of the renormalization group, the 
various quantities in the above equation scale as [l9[ 
r(l*) = re r / u , u{l*) = u* , and <f>{l*) = ^ d ~ 2+ ^ v 72 . 
To account for the finite thickness of the wetting films, 
we choose a rescaling factor e l = L, at which point 
the Landau- Ginzburg parameters are r(l*) = rL x l v and 
u(l*) w u* . Ignoring fluctuations at scales larger than 
(and equal to) L is then equivalent to carrying out the 
previous mean-field calculations with the above rescaled 
parameters (and with L' = 1). Any free energy density 
(and hence effective force) in this approximation is thus 
obtained from its mean-field form as 

In particular, from Eq. (|I0[). the effective force for — ir 2 < 
y < becomes 
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(13) 



The above mean-field approach with renormalized pa- 
rameters, also fixes the vertical scale in Fig. [21 remov- 
ing the unknown value of u. This is because, in units 
where fc^T c = 1, the fourth order coupling constant u* 
has a specific value within renormalization group, and 
is computed to high precision in Ref. [20(. Making use 
of the values provided there [2 1], one finds u* — 0.88. 
Using Eq. (|13[) . the minimum of the scaling function at 
= — 7r 2 is estimated as 
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-6.92, 



(14) 



which is roughly five times larger than the experimentally 
measured value! 

The mean-field calculation suffers from a number of 
other shortcomings. We already noted that the force is 
asymptotically zero for y — > —00 due to the neglect of 
phase and capillary fluctuations It also yields a zero 
critical Casimir force for all y > 0, while there is a small 
finite value to this quantity in experiments due to order 
parameter fluctuations, as calculated at one and two loop 
level in Ref. 0. As a first correction to the mean- field 
result the Casimir force can be estimated by including 
quadratic (i.e. one-loop) fluctuations around the saddle 
point. At exactly y = 0, the fluctuations are mass-less 



4 



and lead (in three dimensions) to a Casimir amplitude 
of @,[22| 



*(0) 



1 f°° 

— / dq q 2 (cothg - 1) = -0.0956566. (15) 
2tt Jo 



We note that the ratio of the extremum value of tt 2 ) 
(a mean field result) to its one-loop magnitude at y = 
is around 70. This is to be contrasted with ratios from 
experimental data that range from 20 to approximately 
35 (the variation of possible ratios arises from the differ- 
ent values of $(0) that can be inferred from experimental 
data) . 

Another feature of the mean-field result is the discon- 
tinuous slope of the scaling function at its minimum in 
Fig. [21 This is a consequence of the onset, in the mean 
field approximation, of an actual transition accompanied 
by long-range order in the film at y = — tt 2 . Calculations 
involving lower dimensional models for the "bulk" and 
"film" (e.g. two dimensional Ising bulk and a one dimen- 
sional Ising film or one dimensional bulk and zero dimen- 
sional film) yield behavior that is qualitatively similar to 
the mean field picture notwithstanding that a transition 
within the film is precluded 22] . In these cases — as may 
well be the case in the experimental setups of [l[ and 
0| — the free energy of the film begins to decrease due to 
the formation of local pockets of order. In the absence 
of a genuine transition in the film there is no singularity 
at the minimum. In the actual films, the onset of super- 
fluidity (which should occur at a different temperature 
altogether) will be of the two-dimensional XY universal- 
ity class 0, [24)]. These details are certainly not taken 
into account in our treatment and, arguably, are not rel- 
evant for the current level of experimental resolution. 

Indeed, let us reemphasize that the essence of our 
derivation is thermodynamic: outside a narrow "critical 
window" thermodynamic signatures of ordering (whether 
ultimately short or long-ranged) should be accurately de- 
scribed by a Landau theory. A phenomenological calcula- 
tion of the thinning effect may be performed on this basis 
using standard mean field approximation, The predicted 
thinning of the film is consistent with the experimental 
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results s 

In summary, we find that a mean field calculation 
yields results for the scaling form of the Casimir force 
in 4 He films that are in qualitative agreement with re- 
cent experimental observations. While discrepancies re- 
main to be resolved, we are confident that this approach 
captures the important thermodynamic signatures, and 
therefore the effective thinning force in 4 Hc films at the 
onset of bulk superfluidity. Work on an improved mean 
field approaches and the proper evaluation of the ef- 
fects of fluctuations in all temperature regimes is ongo- 
ing. An understanding the underlying principles behind 
the thinning of Helium films is expected to have impor- 
tant implications in the analysis of other wetting experi- 
ments [25j, |26j. 
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